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ABSTRACT 

We use the Herschel SPIRE color-color diagram to study the spectral energy distribution (SED) and the redshift estimation of high-z 
galaxies. We compiled a sample of 57 galaxies with spectroscopically confirmed redshifts and SPIRE detections in all three hands 
at z = 2.5 — 6.4, and compared their average SPIRE colors with SED templates from local and high-z libraries. We find that local 
SEDs are inconsistent with high-z observations. The local calibrations of the parameters need to be adjusted to describe the average 
colors of high-z galaxies. For high-z libraries, the templates with an evolution from z=0 to 3 can well describe the average colors of 
the observations at high redshift. Using these templates, we defined color cuts to divide the SPIRE color-color diagram into different 
regions with different mean redshifts. We tested this method and two other color cut methods using a large sample of 783 Herschel- 
selected galaxies, and find that although these methods can separate the sample into populations with different mean redshifts, the 
dispersion of redshifts in each population is considerably large. Additional information is needed for better sampling. 
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Introduction 


|2014| l. Consequently, it is certainly useful to directly estimate 
redshifts for large samples of galaxies only detected in IR. 


^ The Herschel Space Observatory, with its 3.5 meter mirror and 

_the two PACS ( |Poglitsch et al. |20101 l and SPIRE ( [Griffin et al.| 

|2010| l instruments, covers the electromagnetic spectrum from 55 
m to 672 fim. Herschel provides us with an unprecedented chance 
to study the infrared (IR) emission of galaxies. The SPIRE in- 
strument is able to detect galaxies up to quite high redshift (z ^ 3 
to z ^ 7), where the 250, 350, and 500 pm bands of SPIRE corre- 
^ sponds to the peak of the IR spectral energy distribution (SED) of 
star-forming galaxies. The combination of the three band fluxes 
T—( provide constraints on the shape of IR SEDs of these high-z ob- 
ILJ jects and on physical parameters like the dust temperature or the 


dust mass (e.g., Amblard et al.|2010| Chapman et aI.|20I0[|Mag 


[nelli et al.|2010 ]Roseboom et al.|2()i2|l. 




The identification of IR sources and the measure of their red¬ 
shift are challenging at z > 2.5. The obscuration of these high- 
z IR luminous galaxies makes their optical emission very faint 
and hard to detect. The large beam size of SPIRE also makes 
the cross-matching with optical sources very difficult. As a con¬ 
sequence, acquiring spectroscopic or photometric redshifts from 
optical bands is not possible for most of these sources. Submil¬ 
limeter spectroscopic redshifts can be used for these high-z IR 
sources: the redshifts of most spectroscopically confirmed high- 
z sources are measured in this way (e.g., Riechers et al. 2013[ 
Dowell et al. 2014|l. However, the method is not yet efficient 


The IR data from Herschel alone are not sufficient to provide 
reliable photometric redshifts for individual galaxies, but they 
can be useful in providing redshift distributions of large samples 
and candidates for follow-up spectroscopic observations. Several 
works have been dedicated to presenting methods to select high- 
z candidates using the SPIRE data (e.g., [Amblard et alT 2010[ 


[Pope & Chary|20l0{ [Roseboom et a\.\20\2\ . Based on the red¬ 

shift distributions of IR galaxies on a SPIRE color-color dia¬ 
gram, a color cut is applied to select high-z and low-z galaxy 
populations (e. g., [Amblard et'af [2010[ Pope & Chary[ 2010[ 


[Riechers et al.[[20i3[ [Dowell et al.[[2014 i. All these works are 

based on modified blackbody (MBB) models to model the IR 
emission. It is necessary to reinvestigate this topic with updated 
libraries of more sophisticated IR SED templates, which are now 
available using Herschel observations. 

Empirical template of IR SEDs (8 — lOOOqm) are widely 
used to derive the total IR luminosity, the SER, and the dust prop¬ 
erties of galaxies. Prior to Herschel, local templates were built 
based on IRAS, ISO, SCUBA, or Spitzer (e.g.. Chary & EIbaz[ 


[20m][DaIe & HeIou[2002[[Draine & Li[2007l[Rieke et al.[2009| l. 

Adding Herschel data, recent studies have constructed new IR 
templates for low- and high-redshift objects ([Elbaz et al.|2011[ 


enough to obtain redshifts for large galaxy samples ([Casey et al.[ 


Smith et al.|2012b][Magdis et al.[201^ Ciesla et al.|2014[ l. How- 

ever, there are still no IR SED templates for galaxies at z > 3. 
Numerous studies still use local templates or MBB models when 
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deriving the properties of very distant galaxies (e.g., Cox et al.| 

201 1| Conley et al. 

2011 Combes et al.[2012| Huang et al. 2014 

Robson et al.|2014 



To obtain a clearer and up-to-date view of where we stand in 
this post-Herschel era, we compare the SED templates from dif¬ 
ferent libraries with a sample of high-z galaxies that have spec¬ 
troscopic redshifts and reliable SPIRE flux measurements. We 
attempt to find the SED templates that can best describe the av¬ 
erage SPIRE colors of high-z galaxies. Through a comparison 
with the local calibrations, we can test the evolution of the galaxy 
populations. In addition, we examine whether different redshift 
populations could be separated by using a SPIRE color-color di¬ 
agram. 

The paper is organized as follows. We first introduce the 
samples and the SED libraries used in this study in Sections 
and[^ We compare SED templates with observed high-z data on 
a color-color diagram at different redshifts in Section]^ In Sec- 
tion|^we discuss the possible influence of the detection limits of 
SPIRE observations and the ability of applying SPIRE color cuts 
to select high-z galaxies. The main conclusions are summarized 
in Section |6] 

Throughout this paper, 5;^ [mJy] means the flux density Sy at 
A /im. The IR luminosity Lir [Lq] is integrated from 8 to 1000 
/rm of the spectrum. The A-cosmology is adopted: f2,„=0.3, 
Q.a=0.1, and Ho=7Q km s^^ Mpc^*. 


2. Sample 

2.1. A spectroscopically confirmed sample of galaxies at 

z>2.5 


We select from the literature a sample of 57 galaxies at z > 2.5. 
We attempt to include all galaxies that are spectroscopically con¬ 
firmed and have detections in all SPIRE bands in this sample, 
excluding those with apparent signatures of active galactic nu¬ 
cleus (AGN). The galaxies and their sources are listed in Table 
[T] Since a considerable number of these galaxies are lensed, the 


The 

tiere- 


magnification factor of each galaxy is also listed in Table 
galaxies in this sample are taken from different works, and t 
fore are inhomogeneous in selection methods as well as redshift 
measurements. We shortly describe each work in the following. 
\WeiB etal. i2013): The SPT sample from WeiB et al. P013| l 


includes 1.4mm selected galaxies in the South Pole Te 
Sunyaev-Zeldovich (SPT-SZ) survey ([Carlstrom et al. 


escope 

2orr|i 


Synchrotron dominated sources classified by Vieira et al. ( |20l6| l 
have been removed from this sample. The spectroscopic red¬ 
shifts of these galaxies are derived from the Atacama Large Mil¬ 
limeter/submillimeter Array (ALMA) observations. Sixteen of 
these sources are at z > 2.5 and have SPIRE detections in all 
three bands. For three out of the 16 galaxies, Hezaveh et al. 
( 2013| l derived the IR luminosities using the MBB model with 
p = 2.0, and their lens magnification factors from lensing mod¬ 
els. There is no measurement of the Ljr or lens magnification 
factor available in the literature for the left sources. We assume 
no lens magnification for these galaxies in Table [1] however, 
these galaxies are likely to be lensed. 


\MagneUi et al.\ (201 2P : [Magnelli et al.| ( |2012[ ) collected a 
sample of submillimeter galaxies with spectroscopic redshift and 
Herschel observations. Sources with no PAH features are consid¬ 
ered to be AGN and were excluded from this sample. We select 
11 galaxies in their sample with z > 2.5 and detections in three 
SPIRE bands. In their work, the Lir is calculated from multitem¬ 
perature SED fitting. 
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[Casey et aJ.| (2012^ : Six out of the 36 Herschel sources of 
the sample Casey et aL] ( 2012| l are at z > 2.5 and have no sign 
of AGNs according to their optical spectra. The redshift of these 
galaxies was spectroscopically confirmed with the Keck 1 Low 
Resolution Imaging Spectrometer (LRIS) and the Keck II DEep 
Imaging Multi-Object Spectrograph (DEIMOS). The Lir are ob¬ 
tained fitting a single MBB with j3 = 1.5 and a mid-infrared 
(MIR) power law with a slope of a = 2.0 to the IR data. 


Saintonge et al. 


(201 3p : [Saintonge et al.j ([2013|l present a 
sample of lensed galaxies at z = 1.4 - 3.1 with HerscheUSPlRE 
data. These galaxies have intrinsic properties similar to ultra¬ 
violet (UV)-selected lyman break galaxies (LBGs). A compari¬ 
son of these lensed galaxies with unlensed galaxies suggests that 
they are representative of star-forming galaxies at their redshifts. 
The [NII]/Ha ratios of these galaxies are in the range expected 
from star-forming galaxies, indicating that they do not harbor 
AGNs. Six of these galaxies are at z > 2.5 and have detection in 
each SPIRE band. We remove one galaxy (J1527) from our sam¬ 
ple because its SPIRE photometry has very large uncertanties 
(g{Sv)/Sv ~ 10) and is not consistent with its PACS measure¬ 
ment (see their Figure 5). The IR luminosities of these galaxies 
given by [Saintonge et al. ( 2013| l are calculated from the SED 
fitting using [Draine & Li ( 2007| l models. 

Dowell et al. (201 41: [Dowell et ^ (2014[) selected a sample 


of ‘500 pm risers’, which are high-z galaxies, as confirmed by 
the follow-up observations with CARMA. The IR luminosities 
for these galaxies are calculated by fitting MBB models (leaving 
j3 as a free parameter) to the SPIRE photometry and assuming 
no lensing magnification. We assume the AGN contribution is 
negligible. 


Roseboom et aJ.[(2Qj2p : Four of the submillimeter sources 
of [Roseboom et a~ (2012| l have detections in all three SPIRE 
bands. They used an MBB model with j3 = 1.8 and constant 
optical depth to fit the SED of the sources and derive their IR 
luminosity, assuming AGN contribution is negligible. 


Negrello et al. ^201 Op: The two galaxies of Negrello et al. 


( |2010| l, SDP81 and SDP130, are selected using a flux density 
cut at 500 pm (Ssoo > 100 mJy). They are strongly lensed with 


a magnification factor of 19 and 6, respectively. Their total in- 

trinsic Lir are derived by[Negrello et al. 

(2010[l from multi wave- 

length SED fitting using models from 

da Cunha et al. (2008|l. 


The MIR emission of the two galaxies suggests little AGN con¬ 
tribution (Negrello et al.|20r4|. 


[SmoJcic etaJl(2Q7^ : AzTECl and Vd-17871 have z > 2.5. 
The redshift of AzTECl is confirmed to be 4.34 by SMA obser¬ 
vation (Yun et al. in prep.). The redshift of Vd-17871 is obtained 
by DEIMOS/VLT follow-up observation. The Lir given in this 
work is calculated by fitting the SEDs using Draine & Li ( 2007| l 
models. The widths of Lya lines of the two galaxies rule out the 
possibility that they are powerful AGNs. 


Sldias et aJ.|(20i4|): We select A68/nn4 from the work of 


Sklias et al. (2014i, who present Herschel observations for sev¬ 


eral lensed star-forming galaxies with well-determined spectro¬ 
scopic redshifts at z = 1.5 — 3. The object A68/nn4 is at z = 3.19. 


The Lir in this work is calculated from SED fitting using Chary 
[& Elbaz[ (2001;CE01) templates. 

SPT0538-50: SPT0538-50 is a star-forming galaxy first dis¬ 
covered by the South Pole Telescope. The redshift is measured 
by Aravena et al. (2013| l from CO observation. Bothwell et al. 
(2013|l present the Herschel observation for this galaxy and caT 
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Yuan et al. (2014): The SED of galaxies at z > 2.5 


Table 1. High-z galaxies at z > 3 from literature. 


Reference 

Objects 

z 

5250 [mJy] 

^350 [mJy] 

^500 [mJy] 

log(LiR[L0j)" 


Weiss+2013 

SPT0103-45 

3.0917 

121.0+15.0 

210.0+23.0 

222.0+24.0 

13.84 

1.0 

Weiss+2013 

SPT0345-47 

4.2958 

242.0+25.0 

279.0+29.0 

215.0+23.0 

14.3 

1.0 

Weiss+2013 

SPT0346-52 

5.6559 

136.0+16.0 

202.0+22.0 

194.0+21.0 

14.36 

5.4 

Weiss+2013 

SPT0418-47 

4.2248 

115.0+14.0 

189.0+20.0 

187.0+20.0 

14.08 

21.0 

Weiss+2013 

SPT0441-46 

4.4771 

62.0+10.0 

98.0+12.0 

105.0+13.0 

13.87 

1.0 

Weiss+2013 

SPT0459-59 

5.7993 

35.0+10.0 

54.0+10.0 

61.0+11.0 

13.83 

1.0 

Weiss+2013 

SPT0529-54 

3.3689 

74.0+13.0 

137.0+17.0 

162.0+19.0 

13.75 

9.4 

Weiss+2013 

SPT0532-50 

3.3988 

214.0+23.0 

269.0+28.0 

256.0+27.0 

14.08 

1.0 

Weiss+2013 

SPT2103-60 

4.4357 

43.0+10.0 

72.0+11.0 

108.0+15.0 

13.81 

1.0 

Weiss+2013 

SPT2146-55 

4.5672 

58.0+12.0 

79.0+14.0 

82.0+14.0 

13.8 

1.0 

Weiss+2013 

SPT2147-50 

3.7602 

73.0+12.0 

114.0+14.0 

116.0+15.0 

13.77 

1.0 

Weiss+2013 

SPTOl 13-46 

4.2328 

22+8 

54+10 

82+11 

13.63 

1.0 

Weiss+2013 

SPT0125-47 

2.51480 

785+79 

722+73 

488+50 

14.37 

1.0 

Weiss+2013 

SPT0243-49 

5.699 

18+8 

26+8 

59+11 

14.06 

1.0 

Weiss+2013 

SPT2132-58 

4.7677 

55+11 

75+12 

78+12 

13.81 

1.0 

Weiss+2013 

SPT2134-50 

2.7799 

346+36 

339+35 

257+28 

14.09 

1.0 

Magnelli+2012 

LESSOll 

2.679 

15.2+2.6 

19.2+3.4 

17.9+4.2 

12.73 

1.0 

Magnelli+2012 

LOCK850.14 

2.611 

19.4+3.8 

24.9+5.7 

30.6+8.3 

12.82 

1.0 

Magnelli+2012 

AzLOCK.lO 

2.56 

12.4+3.9 

34.5+7.5 

35.0+8.7 

12.8 

1.0 

Magnelli+2012 

SMMJ14011-b0252 

2.565 

61.7+6.0 

63.1+3.8 

48.5+4.0 

13.27 

3.0 

Magnelli+2012 

SMMJ14009-b0252 

2.934 

66.4+6.0 

65.8+3.8 

53.7+4.2 

13.41 

1.5 

Magnelli+2012 

SMMJ00266-bl708 

2.73 

52.8+7.2 

61.3+4.5 

44.6+4.7 

13.28 

2.4 

Magnelli+2012 

SMMJ04542-0301 

2.911 

76.0+15.1 

94.4+8.8 

84.1+11.5 

13.5 

50.0 

Magnelli+2012 

SMMJ16354-I-66114 

3.188 

20.5+6.0 

25.7+3.0 

17.4+3.3 

13.03 

1.7 

Magnelli+2012 

SMMJ16355-I-66122B 

2.516 

78.1+6.0 

70.2+7.2 

74.7+3.5 

13.34 

22.0 

Magnelli+2012 

SMMJ16355-b66123A 

2.516 

53.7+6.0 

47.1+3.1 

26.9+3.5 

13.23 

14.0 

Magnelli+2012 

SMMJ02399-0136 

2.81 

65.7+6.1 

72.0+3.9 

62.8+4.6 

13.38 

2.5 

Casey+2012 

X24 J033136.96-275510.9 

3.145 

15.5+3.8 

20.8+3.8 

24.0+4.8 

12.91 

1.0 

Casey+2012 

X1.4 J033151.94-275326.9 

2.938 

12.9+3.9 

21.6+3.7 

24.2+4.2 

12.83 

1.0 

Casey+2012 

X24 J095916.08-I-021215.3 

4.454 

25.8+2.2 

24.1+2.9 

14.7+3.2 

13.3 

1.0 

Casey+2012 

X1.4J100111.52+022841.3 

3.975 

24.5+2.2 

32.9+4.3 

22.8+6.8 

13.27 

1.0 

Casey+2012 

X24 J100150.16+024017.2 

2.883 

13.8+2.2 

12.2+2.9 

8.9+4.0 

12.71 

1.0 

Casey+2012 

X1.4J100008.64+022043.1 

2.888 

16.5+2.2 

12.3+2.9 

9.4+3.6 

12.79 

1.0 

Saintonge+2013 

8:00arc 

2.73 

91.3+3.6 

78.5+4.2 

47.8+4.0 

13.5 

12.3 

Saintonge+2013 

cB58 

2.729 

24.9+2.4 

16.0+3.3 

10.9+3.6 

12.97 

31.8 

Saintonge+2013 

J0712 

2.646 

9.1+2.9 

3.8+3.4 

4.7+3.1 

12.5 

27.5 

Saintonge+2013 

J1226 

2.925 

9.7+5.3 

3.7+3.3 

3.1+5.1 

12.49 

40.0 

Saintonge+2013 

Eye 

3.074 

26.8+7.0 

33.1+8.7 

19.1+10.0 

13.08 

30.0 

Dowell+2014 

LSW20 

3.36 

17.6+0.0 

36.6+0.0 

43.9+0.0 

13.17 

1.0 

Dowell+2014 

ELSl 

4.29 

62.3+0.0 

94.0+0.0 

98.2+0.0 

13.81 

1.0 

Dowell+2014 

ELS5 

4.44 

24.7+0.0 

45.1+0.0 

48.2+0.0 

13.5 

1.0 

Dowell+2014 

LSW102 

5.29 

49.7+0.0 

118.1+0.0 

140.4+0.0 

14.11 

1.0 

Roseboom+2012 

G1 

4.05 

18.6+2.3 

31.7+2.7 

24.3+2.7 

13.23 

1.0 

Roseboom+2012 

G3 

2.98 

46.3+2.0 

46.5+2.3 

32.7+2.3 

13.28 

1.0 

Roseboom+2012 

G19 

2.91 

11.2+2.1 

11.3+2.4 

17.3+2.4 

12.65 

1.0 

Roseboom+2012 

LI 

2.562 

80.1+3.4 

63.2+3.8 

40.6+4.4 

13.4 

1.0 

Negrello+2010 

SDP81 

3.04 

129.0+20.0 

182.0+28.0 

166.0+27.0 

13.79 

19.0 

Negrello+2010 

SDP130 

2.63 

105.0+17.0 

128.0+20.0 

108.0+18.0 

13.54 

6.0 

Smolcic+2015 

AzTECl 

4.3415 

19.5+6.0 

29.8+7.6 

28.8+9.0 

13.31 

1.0 

Smolcic+2015 

Vd-17871 

4.622 

9.5+2.0 

12.3+2.4 

13.0+3.7 

13.01 

1.0 

Sklias+2014^ 

A68/nn4 

3.19 

45.3+1.87 

38.04+2.19 

17.74+1.5 

13.01 

2.3 

Bothwell+2013 

SPT0538-50 

2.8 

326.0+23.0 

396.0+38.0 

325.0+24.0 

14.09 

21.0 

Magdis+2011 

GN20 

4.055 

18.6+2.7 

41.3+5.2 

39.7+6.1 

13.34 

1.0 

Cox+2011 

HATLASID141 

4.24 

115.0+19.0 

192.0+30.0 

204.0+32.0 

14.1 

20.0 

Conley+2011* 

HLSW-01 

2.957 

425.0+10.0 

340.0+10.0 

233.0+11.0 

14.22 

11.0 

Combes+2012 

HLSA773 

5.24 

85.0+8.0 

168.0+8.0 

203.0+9.0 

14.26 

11.0 

Reichers+2013^ 

HELS3 

6.34 

12.0+2.3 

32.4+2.3 

47.3+2.8 

13.79 

<3.5 


Notes. IR luminosity derived using DH02 template fitting, uncorrected for lensing magnification. Calibration errors and confusion noise are 
not included. 
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log(L|p,(literature)[LJ) 


These high-z galaxies are bright in IR, suggesting that they 
are starbursting galaxies. However, at high-z the definition of 
main-sequence (MS) galaxies and starbursts (SBs) is not clear. 
Several works have shown that the extremely high star forma¬ 
tion rate (SFR) of high-z sources is due to the large amount of 
cold gas, and that SBs should be defined according to the com¬ 
pactness of the star-forming region instead of the SFR ( |Elbaz| 
et al.| 20 lT]l. In the case of our sample galaxies, further informa¬ 


tion is needed to be able to state if they belong to the MS or not. 
Therefore, we do not distinguish MS sources from SBs when we 
compare them with SED libraries in Section ]?!^ 

In our sample we have removed galaxies with apparent AGN 
features, however, we cannot rule out galaxies with less power¬ 
ful or obscured AGNs. Moreover, the AGN identification is not 
uniform due to the inhomogenous data used in each work. As a 
result, some galaxies in our sample may have embedded AGNs. 
In Section 3.4 we show that the AGN contribution is insignificant 
for the wavelength range discussed in this work. 


Fig. 1. Comparison of the observed IR luminosity Ljr derived from 
SED fitting based on the templates of DH02 and the values from litera¬ 
ture. 


culated the IR luminosity using CIGALiQsED fitting. The 


& Helou (2002;DH02) models are used for IR part. 


Dale 


GN20: GN20 is first dete cted by|Pope et al.| ( 2006) 1 with a 
redshift z = 4.05 obtained by Daddi et al.| ( |2009 i. The SPIRE 
fluxes were obtained by |Magdis et al. ( |201 1^, who a lso estimated 
the IR luminosity using the models of Draine & Li (2007|l. 


H-ATLAS ID141: H-ATLAS ID141 is a source lying at 
z = 4.24 with an apparent IR luminosity of ~ (8.5 ± 0.3) x 
LQ obtained from an MBB fitting, under the optically 
thin assumption, and a magnification factor between 10 and 30 
dCox et al.| 20 TT] ). 


HLSW-01: HLSW-01 (HERMES J105751 1-H573027) was 
spectroscopically confirmed at z = 2.96 by Riechers et al. 


(201 l|l. This multiple source is lensed by a group of galaxies 


resulting in a magnification factor of ji — 10.9 ±0.7 (Gavazzi 
[eTarilT^ . The SPIRE fluxes of HLSW-0 1 are from the Her- 
MES SCATv3.1 cat alog ([Smith et ah 2012a| l, and its IR luminos¬ 
ity is calculated by |Conley et al. ( 2()1 1| | using SED fitting with 
an MBB model. 


HLSA773: HLSA773 is a source at z = 5.24 with an IR lu¬ 
minosity of 1.1 X Lf. 

MBB model with j3 

[20T2] |. 


.0 obtained from an optically thin 
2, and jj. estimated at ~ 11 (Combes et al. 


HFLS3: HELS3 satisfies ^soo > ^350 > ^ 250 , has a redshift 
of z = 6.34, and a Ljr of 4.16 x 10*^ Lq estimated from an MBB 
SED fitting ( [Riechers et al.|[2013| l. According to |Cooray et al.| 
( 2014| l, the lensing magnification of this galaxy is less than 3.5. 

The IR luminosities from the literature are derived using dif¬ 
ferent wavelength ranges and different SEDs. To present more 
uniform measurements of IR luminosities, we recalculate Ljr 
by fitting the SEDs using DH02 templates and the three-band 
SPIRE fluxes. The results are shown in Table (the values are 
not corrected for lensing magnification). We compared the re¬ 
sults with the literature values in Eigure [T] and found that our 
results are in agreement with the estimations given in the litera¬ 
ture. 


* Code Investigating GALaxy Emission, http://cigale.lam.fr (Noll 


et al. 2009 1 
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2.2. Comparison samples 


In addition to the high-z sample of identified objects, we use 
other samples for comparison. Eirst, a UV-selected sample at 
z 3 and 4 compiled by |Heinis et al.| ( [2014| l from the Herschel 
Multi-tiered Extragalactic Survey key program (HerMES, |01iver| 
et al.||201^. The redshifts of these galaxies are taken from the 


photometric redshift catalogue built from COSMOS data by |IT] 
|bert et al!] ( |2009| version 2.0), and the mean error of the redshift 
is about 0.4 and 0.7 at z ~ 3 and z ~ 4, respectively. The IR fluxes 
are estimated from HerMES data using stacking analysis. The IR 
luminosity Ljr of the stacked sample is ~ 4 x IO^Lq. We use 
the UV-selected galaxies to compare with the spectroscopically 
confirmed high-z sample in Section [431 

When discussing the methods to select high-z populations 
using the SPIRE color-color diagram (Section |5.2| l, we use 
a sample of 500 pm selected galaxies taken from Rowan- 


Robinson et al. (2014p. These sources are obtained by match 


ing the HerMES SCAT catalog ( [Oliver et al.|2012 


|2014[ l in the SWIRE-Lockman area ([Lonsdale et 


Wang et al. 
.|2003[ ) with 


the SWIRE catalog ( [Rowan-Robinson et al.|2008| l, with the re¬ 
quirement that the 24pm fluxes are brighter than 100/fJy. The 
redshifts are derived using up to 13 photometric bands available 
in the Lockman area. When there are multiple counterparts in the 
SWIRE catalog, the source with the highest 450/4m predicted ac¬ 
cording to the SED is selected. The sample include 1335 SPIRE 
500/im selected sources, 967 of which have SWIRE associa¬ 
tions and 783 sources are with detections in all the three SPIRE 
bands. The redshifts of the sample range from z ~ 0 to z ~ 5. The 
photometric redshift error for this sample is about 4% in (1 ±z) 
( [Rowan-Robinson et al.||2OT4| ). The selection effect and the re¬ 
liability of the redshifts for this sample are discussed in Section 


3. SED tracks on the SPIRE color-color diagram 

3.1. Single-temperature modified blackbody models 


We use an MBB model to derive the dust temperature of the 
galaxies of our sample, as commonly done in the literature (e.g., 
Roseboom et al. 2012[ Casey et al.[|2012[ Huang et al.[|2014} . 


with 


/v = (l-e-"^)Bv(7)(), 


( 1 ) 
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where the optical depth Xx (v/vo)^, j3 is the grain emissiv- 
ity index with values ranging from ~ 1.5 to ^ 2 (Hildebrand 
1983| l, Vo is the frequency at which the radiation becomes opti¬ 
cally thick, and ByiTd) the blackbody spectrum at a dust tem¬ 
perature of To. The upper panel of Figure shows the SPIRE 
color-color diagram of the MBB model with dust temperatures 
at 30, 40, and 60 K. Here, Vo is assumed to be c/lOO/fm and 
j3 =1.8. We follow the results of Roseboom et al.| ( |2012| l, who 
found that a T — 40K MBB model with Vq = c/lOO/im and 
j3 = 1.8 can adequately describe the color of their sample of 
galaxies with 2 < z < 5. The upper panel of Figurej^shows that, 
at the same redshift, cold dust galaxies are above and more to 
the left than hot dust galaxies, i.e., at higher 5500/^350 and lower 
Sxso/Ssso compared to hot dust galaxies. Therefore, we refer to 
the upper-left part of the SPIRE color-color diagram as the ‘cold 
region’, while the lower-right part of the diagram is the ‘hot re¬ 
gion’ according to the MBB models. The lower panel of Figure]^ 
shows that taking different value of j3 can cause systematic dif¬ 
ferences when estimating the dust temperatures. However, this 
does not affect our conclusions since we are only interested in 
the relative dust temperatures. In the following part, we assume 
Vo = c/100/rm and )3 =1.8 for MBB models. 


3.2. SED libraries 


The SED libraries used in this work are CEOl, DH02,|Elbaz et al.| 
(20 11; El l),|Magdis et al.| (2012; M12),|BertaeCT(2013; B13), 
and |Ciesla et al.| (2014; C14), shown in FigureMost of these 
libraries are based on data at z < 3. Currently, no library exists 
for z > 3 objects. Therefore, we need to extrapolate these tem¬ 
plates to higher redshift to compare with the data. According to 
M12, at z > 3 the SED of a MS galaxy stops evolving because of 
the flattening of the evolution of the specific star formation rate 
(sSER). Therefore, we assume no evolution after z > 3 for these 
libraries. For z < 3, only the M12 MS library takes the evolu¬ 
tion of IR SEDs with redshifts into account, so we set the tem¬ 
plates of M12 as the reference and compare all other templates to 
them. The color-color diagrams for the templates selected from 
these libraries are shown in Eigure|^ In the following, we intro¬ 
duce the libraries and discuss the position of the templates on the 
SPIRE color-color diagram. 


3.2.1. Pre-Hersche/templates 


Pm-Herschel libraries, such as CEOl and DH02, are commonly 


used in the studies of high-z galaxies (e.g., Magdis et al. 

2010 [ 

Elbaz et al.|201 lUBerta et al.|2013[|Riechers et al.|2013[| 

leinisi 

et al. 

2014 

1 . The templates of CEOl are built to reproduce the 
and SCUBA observations, using the dust models of 

ISO,] 

[RAS 


Silva et ar| ( |1998| l and ISO MIR spectra. These templates are 
parameterized by the total IR luminosity, ranging from Lir ~ 
10^ to ~ 10*^ Lq. The templates of DH02 are based on IRAS, 
ISO, SIRTF, and SCUBA observations. The local SEDs based 
on the models of |Desert et al.| ( |T990| are combined assuming the 
distribution of dust mass as r/Mj 0 = dU, where U is the 
intensity of the radiation field and the exponent a ranges from 
1 to 2.5 for star-forming galaxies. Galaxies with small a values 
are more actively forming stars. 


3.2.2. Post-Hersche/ templates 

Based on GOODS-Herschel data at 0 < z < 2.5 and local 
data of AKARI, ISO, and IRAS, |Elbaz eTakl ( |20TT] ) found that 


MocJifiecj Blackbody (p=1.8) 



Modified Blackbody (T=40K) 



S25C/S350 


Fig. 2. The SPIRE color-color relation for MBB models. Upper panel: 
MBB models with j3 = 1.8, Vq = c/lOO/im and r(]=30K (Blue solid 
line), 40 K (green dashed line), and 60 K (red dotted line). Lower panel: 
MBB models with T =40 K and jS = 2.0 (Blue solid line), 1.8 (green 
dashed line), and 1.5 (red dotted line). The numbers on the plot indicate 
the redshifts. 


most ULIRGs at high-z are actually forming stars in a “quies¬ 
cent” mode similar to local MS galaxies, and their high SFR 
are due to the large amount of molecular gas. Therefore, they 
suggest universal SED templates for MS galaxies and for SB 
galaxies without assuming any evolution. The templates were 
built by fitting their sample (0 < z < 2.5), using a “diffuse ISM” 
component and a “star-forming region” component. The models 
of the two components are given by |Galametz et al.| ( |2009) l and 
|Galliano et al.|(|2011|). Th e galactic dust properties are adopted 
from [Zubko et al. 1^004) 1. The distribution of dust mass heated 
by different starlight intensities U is assumed to be a power law; 
dMi/dU oc f/-“. 


M12 derive the MS SEDs by fitting the Draine & Li ( |2007| l 
model to the data of a well-defined set of normal MS galaxies, 
including nine normal star-forming galaxies at z ~ 0.5 and ^1.5, 
three high-z SMGs, and a stacked sample of > 4000 MS galaxies 
at 0.5 < z < 2.5, all of which have spectroscopic redshifts, rich 
rest-frame UV to MIR photometry and Herschel detection. They 
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low-z 


high-z 




Fig. 3. SED templates studied here. Left panel: low-z templates. In navy: CEOl templates with Ljr ~ 10‘^, 10*®, 10**, 10*^, 10*® Lq. In yellow: 
DH02 templates with a = 1.0, 1.5, 2.0, 2.5. In magenta: C14 mean template. Right panel: high-z templates. In red: Ell templates. The solid line 
denotes the MS template and the dashed line denote the SB template. In black: M12 templates. The solid line denotes the MS template and the 
dashed line denotes the SB template. In green: four SF templates selected from B13. 


find that SEDs of MS galaxies can be characterized by the dust- 
weighted mean starlight intensity scale factor <U >, which is 
proportional to Lir/Mj, and is mildly evolving with cosmic time 
as <U (1 -l-z)' They build eight SED templates (U1 to 
U8) for MS galaxies at eight redshift bins (Table |^. In addition, 
they suggest that SB galaxies have a universal SED, and use the 
ht of GN20 (z ~ 4) as the template for SB galaxies. 


Table 2. Redshifts range for M12 models 


MS templates 

z 

U1 

0 - 0.025 

U2 

0.05 - 0.275 

U3 

0.3-0.625 

U4 

0.65 - 0.975 

U5 

1.0-1.30 

U6 

1.325- 1.725 

U7 

1.75-2.25 

U8 

2.27 -3.0 


|Berta et ah ( |2013| l built 32 SED templates from PEP and 
HerMES data. Each template is the median SED of a class of 
galaxies defined from the multiwavelength color space. Eight of 
these SED templates are for rare objects (outliers of the main 
classes defined from the color space), including torus, AGNs or 
even elliptical galaxies. The SEDs are fitted using MAGPHYS 
code ( |da Cunha et al. 12008) 1. 

Other post-Herschel templates are derived from local obser¬ 
vations (e.g., [Smith et al.||2012bl C14). We use the template 
of C14 to compare with the high-z libraries. Using Draine & 
|Li] ( |2007] i models, C14 fit the UV to sub millimeter SEDs of a 
sample of 322 nearby galaxies, the Herschel Reference Survey 
(HRS, |Boselli et al.|2010| l, which is volume limited and selected 
in K-band. The fitted SEDs are averaged to provide the mean 
SED template. 
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3.3. SPIRE color-color diagram based on different SED 
libraries 

Eigurej^ shows the SPIRE color-color diagram predicted by dif¬ 
ferent SED libraries (see more color-color diagrams in the Ap¬ 
pendix). The M12 MS templates are plotted in each panel for 
comparison. At z < 3, the U1 to U8 templates are used for the 
redshift bins, which they are assigned to by M12. Beyond z = 3, 
the U8 template is used to derive the colors considering the flat¬ 
tening of the evolution of the sSER. 

Panels (a) and (b) of Pigure|^present the SPIRE color-color 
diagram based on the pre-Herschel templates. Panel (a) shows 
that the template of CEOl with Ljr = 1O**L0 produces a SPIRE 
color-color relation similar to the MS templates of M12. Panel 
(b) shows that in the DH02 library, the template with a = 2.0 
agrees well with the MS templates of M12 in the local universe, 
whereas as z increases to 2, a must reduce to ~ 1.5 to be consis¬ 
tent with the colors reproduced by the M12 MS templates. Be¬ 
yond z = 2, a template of DH02 with a value between 1.5 and 
2.0 seems to be in agreement with the MS template of Ml 2. Pan¬ 
els (a) and (b) of Pigure|^also show that the templates of CEOl 
and DH02 do not produce 5500/^350 values as large as predicted 
by the MBB models. 

Panel (c) shows the comparison between the templates of 
Ell with the MS templates of M12 on color-color diagram. Both 
the MS and SB templates of El 1 predict colder colors than M12 
galaxies at z < 6. The difference is due to the difference of the 
Rayleigh-Jeans side between Ell and M12 SEDs (Eigure|^. 

Panel (d) compares the M12 MS template to their SB tem¬ 
plate, i.e., the SED of GN 20. It appears that at every redshift the 
SB template is located in hotter regions of the color-color dia¬ 
gram compared to the MS templates. This is not surprising since 
the star-forming activity in SB galaxies is more intense than in 
MS galaxies and can thus heat the dust more efficiently. 
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(a) CE01 (L,r~10^^[LJ 



S 250 /S 350 


(b) DH02 



S 250 /S 350 


(c) E11 



(d) Ml2 



S250/S350 


S250/S350 


(e) B13 (SF_glx_2) 



^ 250^^350 


(f)C14 



S250/S350 


Fig. 4. SPIRE color ratios based on different SED libraries. The numbers on the plot indicate the redshifts. The SPIRE color ratios based on the 
M12 templates are plotted in each panel as a reference (red line with squares). Panels (a): template of CEOl with Ljr ~ 10*'Lq; (b): template of 
DH02 with a = 2.0 (black) and a = 1.5 (blue); (c) templates of Ell for MS galaxies (black) and SBs (blue); (d) templates of M12 for MS galaxies 
(red) and SBs (black); (e) SF_glx_2 template of B13, which is derived from star-forming galaxies in COSMOS field; (f) template of C14. 
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We select four star-forming galaxy templates from B13. We 
found that, on a color-color diagram, the template SF_glx_2 (de¬ 
rived from star-forming galaxies of the COSMOS field) repro¬ 
duces SPIRE colors in agreement with M12 MS templates at z > 
4 (Panel [e] of Figure |^. The SF_glx_l template (star-forming 
galaxy SED derived from the GOODS-S field) reproduces much 
colder colors than the M12 templates. The MIRex_SF_glx (MIR 
excess star-forming galaxies) and Obs_SE_glx (obscured star¬ 
forming galaxies) templates have hotter SPIRE colors than the 
M12 MS templates (see Eigure[A3|). 


Panel (f) shows that at z ~ 0 the C14 template agrees well 
with the M12 MS template (Ul) on the SPIRE color-color dia¬ 
gram, whereas at high-z the SPIRE colors of the C14 template is 
much colder than the M12 template. 

Many studies have shown that at A > 100/4m the SED can be 
well approximated by a simple MBB (e.g., Gordon et al.|2010 


Davies et al.||2012]l. Indeed, most of the templates can be well 


parameterized by a single temperature MBB at the SPIRE wave¬ 
length. Therefore, MBB models can be used for simplicity. How¬ 
ever, for models with redshift evolution (e.g., M12) ,or those 
with wide peaks (e.g.. Ell), MBB is not a good approximation. 
Also, when going to higher redshifts (z > 2), the observed SPIRE 
bands are moving to shorter wavelength where the dust emission 
becomes more complicated so that the MBB model is over sim¬ 
plified. 


3.4. Contribution from AGNs 


Some libraries include contributions from AGNs, however, sev¬ 
eral works showed that the AGN contribution is negligible 
(Elbaz et al. [2010 201 1[ Casey et al. 20141 at long wave¬ 
lengths (longer than ^30/rm). Furthermore, Hatziminaoglou 
|et d.| ( |2UT0l l showed that the SPIRE colors of AGNs and non- 
AGN, star-forming galaxies are indistinguishable. We also find 
that the host galaxy dominated AGN template of B13 generates 
SPIRE colors very similar to M12 templates. To include an AGN 
contribution, [Dale et ah ( 2014| l provided a modified version of 
DH02 templates in which the contribution of a quasar of differ¬ 
ent intensity can be added. Comparing the [Dale et ^ P014) l 
template containing an AGN emission with its corresponding 
DH02 template using a = 2.0, we find that when the 5-20 /rm 
AGN fraction is less than 50%, the difference of SPIRE colors 
(5250/^350, 5350 / 5500 ) given by the two templates are less than 
(0.1,0.15). In this paper, since we are working with SPIRE col¬ 
ors where the AGNs hardly affect the results, we do not discuss 
the possible AGN contamination. 


3.5. CMB contribution 

At high redshift the impact of cosmic microwave background 
(CMB) on the observed IR SED can be significant. The CMB 
affects the observed dust emission with two competing pro¬ 
cesses: it heats the ISM in galaxies and produces a consider¬ 
able background against which galaxies are observed ( |da Cunha| 
et al.|2013] l. The CMB effect can complicate the comparison be¬ 
tween data and templates. Eortunately, according to the results 
of |da Cunha et al.| ( |2013| l, for the redshift range considered here 
(z < 6.5) the effect of CMB on dust temperatures is ~ ±3K. Eor 
galaxies with lower intrinsic dust temperatures (e.g., 18K), the 
CMB effect on dust temperatures is a function of redshifts. But 
for galaxies with higher intrinsic dust temperatures (e.g., 40K), 
the effect is almost constant at different redshifts. In Section|4]we 
find that our high-z galaxies have dust temperatures higher than 
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30K. Therefore, the CMB effect on these galaxies is negligible 
for our discussion in this paper. 


4. Results 


4 . 1 . Gaiaxies atz> 2.5 on the SPIRE color-color diagram 


The 5250/5350 versus 5500/5350 values for our high-z sample are 
plotted in Eigure]^ The sample is divided into four redshift bins: 
z ~ 3 (2.5 < z < 3.5), z ~ 4 (3.5 < z < 4.5), z ~ 5 (4.5 < z < 5.5), 
and z 6 (5.5 < z < 6.5). The colors given by MBB models with 
different dust temperatures are also plotted. We assume optically 
thick MBB models, with Vo = c/100/4m and j3 =1.8 (see also 
Section [3.1| l. 

In all panels only three galaxies have 5500/5350 colors larger 
than 1.5: the upper left region of the color-color plot is avoided 
by all the observations. The lack of galaxies with these colors 
may indicate that at high-z most galaxies have intrinsic high dust 
temperatures. Before reaching a conclusion we must check if the 
detection limits of the observations could prevent us from select¬ 
ing cold dust objects: this is discussed in Section [5T| The sam¬ 
ple of Saintonge et al. ( 2013| l (gray squares) seems to have hotter 
dust temperatures than the other galaxies, which may be due to 
the difference between UV- and IR-selected galaxies (Section 


In our high-z sample, the redshifts are measured from differ¬ 
ent wavelength. The reliability of optically measured redshifts 
given in [Casey et al.| ( |2012| l and [Roseboom et al.j ( |2012| l de¬ 
pends on the cross-matching methods. These sources are shown 
as crosses in Figure It seems at z 4 these sources have 
higher dust temperatures comparing to other sources, but remov¬ 
ing these sources does not affect our conclusions. 


4.2. Comparison of SED templates with data atz> 2.5 

To quantify the coherence between templates and data, we fit 
each SED template to all the high-z galaxies and calculate the 

value for each galaxy. The mean value of the x^ (< >) 

used as an indicator of the average quality of the fits with a given 
template. During the fit, we assume a confusion noise of 6 mJy 
for the galaxies with no confusion noise estimation. The results 
are shown in Table The colors corresponding to the templates 
used in this comparison (both local and for high-z objects) are 
also plotted in Eigure]^ 

We briefly discuss the case of the local templates from CEOl, 
DH02, and C14: we do not expect a good agreement since the 
local calibrations are known to be invalid at high-z ( jCasey et~ar 
|2014| |Lutz||2014| ). Indeed, Eigure shows that the fixed tem- 
plate of C14 is too cold to describe the high-z galaxies. The 
CEOl and DH02 templates can be adjusted to fit the emission 
of high-z populations. The template of CEOl corresponding to 
Lir to ^ lO^^L© can reproduce the SPIRE observations of our 
high-z sample. However, the local calibration of these templates 
cannot be used to give correct Ljr measurements. Eor the DH02 
templates, the lowest value of < X^ > is obtained for a = 1.5, 
a value lower than the average one found for galaxies at z < 2 
(a = 2, |Buat et al.|20TTl[20T^ . 

The high-z libraries include the templates of El 1, M12, and 
B13. The MS template of Ell predicts no evolution with red¬ 
shift of the SED shape for MS galaxies. With this template 
< x"^ >= 22.34 indicating the poor quality of the fit. Figure 
shows that the Ell MS template produces colors colder than 
most observations at z > 3, implying that the template is bi¬ 
ased to galaxies with cold dust temperatures. The SB template 
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0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 2.0 


^250/^350 

Fig. 5. Comparison of the high-z observations and the SED templates on the SPIRE color-color diagram in different redshift bins (z ~ 3, 4, 5 and 

6) . The observations data are in gray and the results of the templates are in color. The SED templates plotted in the figure include: 1) MBB models 
with T = 30, 45 and 60K (navy dots); 2) templates of CEOl (blue open diamonds, the sizes of which are increased with the IR luminosity); 3) 
templates of DH02 (red line) with a = 1.0, 1.5, 2.0, and 2.5 marked in triangles with increasing sizes; 4) MS and SB templates of Ell (cyan 
circles and plus signs); 5) MS and SB templates of M12 (magenta squares and asterisks); 6) four selected templates of B13 (yellow open triangles); 

7) the templates of C14 (green crosses). The data plotted in the figure include galaxies from Table[T] The data from ||Saintonge et al.H2013t are 
shown in gray squares, those from |Casey et al.|P012[ l and |Roseboom et aL] ( |2012^ are shown in gray crosses, and the data from other works are 
shown in gray dots. 


of Ell gives a better fit to the observation data, but is still 


colder than most observations and with a high value of < > 

« ^2 >= 13 . 79 ). 
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Table 3. Table of templates and the mean of each template to fit the data at different redshift bins. 


templates 

all 

2.5 <z<3.5 

<X^> 
3.5 <z <4.5 

4.5 <z<5.5 

5.5 < z < 6.5 

local libraries 

CEOl lO'^L© 

29.61 

42.62 

12.53 

7.52 

6.43 

CEOl 10"L© 

7.09 

9.40 

4.24 

1.93 

3.97 

CEOl lO'^L© 

4.78 

5.87 

3.97 

2.62 

1.09 

DH02 a = 1.0 

6.01 

7.04 

5.33 

4.62 

1.32 

DH02 a = 1.5 

4.14 

5.52 

2.78 

0.86 

1.31 

DH02 a = 2.0 

19.77 

28.21 

8.83 

4.82 

5.08 

DH02 a = 2.5 

39.67 

55.84 

18.88 

12.74 

8.69 

C14 

32.56 

44.22 

17.28 

13.87 

10.33 

high-z libraries 

Films 

22.34 

33.11 

7.99 

3.80 

4.11 

Ellsb 

13.79 

20.37 

5.11 

1.86 

3.05 

Ml 2ms 

6.72 

8.94 

3.94 

1.63 

3.94 

M12sb 

4.85 

6.50 

2.95 

0.64 

2.81 

B13SF 1 

32.56 

40.59 

21.70 

21.57 

16.06 

B13 SF_2 

9.78 

13.59 

4.61 

2.13 

4.98 

B13 0bs_SF 

6.87 

7.80 

6.37 

6.17 

1.59 

B13MIRex SF 

9.03 

13.35 

3.54 

1.00 

1.66 


Assuming a mild evolution of the MS SEDs, the M12 tem¬ 
plates show a better consistency with our high-z sample (< 

>= 6.72). Here we use the U8 template, which is assigned 
to z ~ 3 by M12, assuming no further evolution for redshifts > 
3. Figure l^shows that the U8 template is able to reproduce the 
average observed colors at z ^ 4, 5, and 6, in agreement with 
M12’s assumption that the SEDs evolve very little at z > 3 be¬ 
cause of the flattening of the evolution of the sSFR. The SB tem¬ 
plate of M12 gives smaller values of < X^> (4.85) than the MS 
template. M12 used G20 to derive their SB template. Removing 
this galaxy from the high-z sample, the < X^> becomes 7.35 for 
the MS template, and 5.22 for the SB template. The SB template 
still fits slightly better than the MS template, implying most of 
our high-z galaxies more likely belong to SB population. 

The < X^ > values obtained with the B13 models are inter¬ 
mediate. The colors of the B13 template for star-forming galax¬ 
ies in GOODS held (sf_glx_l) extend to colder regions than ob¬ 
servations. The bestfit template in the B13 library is the template 
for obscured star-forming galaxies (Obs_SF_glx), which gives 
the <^2 > of 6.87. 

We also calculate the < different redshift bins. The 

result for each template is shown in Table and Figure The 
< X^ > for sach template at z 6 becomes similar, which means 
that the goodness of the fit becomes less dependent on the tem¬ 
plates because of the large uncertainties of the data. The large 
uncertainties of the data at z ~ 6 causes the < X^ > value of the 
fitting at z ~ 6 to be smaller than at z 3. The the < X^> value 
are smaller at z ~ 6 also because there are only four sources at 
z ~ 6 and thus less scatter between the sources. 

CEOl template with Lir to ^ 10"L©, DH02 template with 
a = 1.5 and M12 template fit the data well at each redshift bin. 
It is possible that galaxies with these kinds of templates are ma¬ 
jor populations at high-z, but the better fitting can also be due to 
the data not being sufficient to constrain the SED. The local tem¬ 
plates CEOl and DH02 give similar or better < X^ > than high-z 
templates, probably because they have flexible parameters {a in 
DH02and Lir in CEOl). 

Since the SPIRE beam is larger at longer wavelengths, the 
blending affects 500/4m fluxes more than 350/4m and 250/4m 
fluxes. Therefore, the observed SEDs may be overestimated at 
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longer wavelength. If the effect is non-negligible, the intrinsic 
colors of the galaxies is bluer than those plotted in Figure]^ The 
templates of CEOl with Ljr > IO^'L©, DH02 with a < 1.5, the 
SB template of M12, and the Obs_SF_glx template of B13 are 
still preferential. 


4.3. UV-selected galaxies on the SPIRE color-color diagram 


Saintonge et al.| ( |2013| l have shown that their galaxies are 
ilar to UV-selected LBGs, with a lower Lir comparing with 
the IR/submillimeter-selected galaxies. Figure shows that the 


galaxies from Saintonge et al. ( 2013| l are located in hotter re¬ 
gions than the other IR/submillimeter-selected galaxies (based 
on the MBB models), probably indicating hotter dust tempera¬ 
tures for these UV-selected galaxies. 

From a stacking analysis in the COSMOS held, |Heinis et al. 


( |2014| l measured average IR fluxes for samples of UV-selected 
galaxies at z ~ 3 (20849 galaxies) and ^ 4 (6794 galaxies) (see 
also Section 2.2 for the details of this sample). Comparing these 
samples with our high-z sample and the sample of |Saintonge| 


et al. ( |2013 > in Figure [7} w e And that, as with the sample of 
galaxies in [Saintonge et ah (2013|l, the SPIRE colors of these 


UV-selected samples are slightly hotter than the average colors 
of our high-z sample of IR/submillimeter-selected galaxies. 

A high dust temperature in a UV-bright galaxy can be caused 
by a strong UV radiation intensity, which can heat the dust grains 
in the star-forming regions to high temperatures. The high dust 
temperatures and low-IR luminosities of the UV-selected sam¬ 
ple as compared to those of the IR/submillimeter-selected galax- 


low {Md oc LirL^ ^+^) ( 

Casey et al. 

2014 

1 . In IR/submillimeter- 

selected galaxies, the radiation hel 

d may evolve less than es- 


timated by previous studies, which use UV-selected galaxies 
to probe their evolution (e.g., Saintonge et al. 2013| l. A de¬ 
tailed comparison of high-z UV-selected and IR/submillimeter- 
selected galaxies is beyond the scope of this paper and will be 
addressed in a future work. 

In summary, from the comparison between the data and tem¬ 
plates we And the following; 1) most galaxies of our sample (ex- 
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CE1 CE2 CE3 DH1 DH2 DH3 DH4 C14 E11ms E11sb M12ms M12sb B13a B13b B13c B13d 

Model 


Fig. 6. The mean for templates listed in|^at z ~ 3 (blue line with asterisks), z ~ 4 (green line with diamonds), z ~ 5 (yellow line with triangles), 
and z ~ 6 (red line with squares). The mean x^ for all the galaxies are plotted in black line with circles. In the figure, we use CEl, CE2, and CE3 to 
represent CEOl templates with 10 ®Lq, IO^'L© , and IO'^Lq, respectively. We use DHl, DH2, DH3, and DH4 to represent DH02 templates with 
a = 1.0, 1.5, 2.0, and 2.5, respectively. We use B13a, B13b, B13c, and B13d to represent the SF_glx_l, SF_glx_2, Obs_SF_glx, MlRex_SF_glx 
templates in the B13 library, respectively. 


cept 3 of them) have ^ 500/5250 < 1.5, inconsistent with MBB 
models with cold dust temperatures or local mean SED templates 
(C14); 2) to fit the observations at high-z, the parameter Lir of 
the CEOl templates should be adjusted to ~ lO'^L©, and the pa¬ 
rameter a of the DH02 template should be adjusted to ~ 1.5; 
3) the Ell MS template corresponds to colder colors than the 
average of the observations; 4) the MS template of M12 can be 
used for high-z galaxies (z > 2.5) to reproduce the average col¬ 
ors given by the observations; and 5) the colors of UV-selected 
galaxies correspond to hotter dust temperatures than those of 
submillimeter/IR-selected galaxies. 

5. Discussion 

5.1. Detection limits 

In Section|4]we showed that there are no observations within the 
upper left regions of the SPIRE color-color diagram, which is in¬ 
consistent with a number of SED libraries. This can be due to the 
intrinsic high dust temperature of high-z sources and/or a selec¬ 
tion effect due to the detection limits. In this section, we check 
whether the detection limits prevent us from observing objects 
in the upper left regions of the SPIRE color-color diagram. 

Eirst, it is possible that the high-z sample is biased against 
galaxies with certain types of SEDs, because a galaxy only 
appears on the color-color diagram if it is detected in the 
three SPIRE bands (color selection effect). Therefore, we check 
whether the detection limits of our high-z sample have an im¬ 
pact on the type of SEDs we can probe. We examined five SEDs 
that are less favored by the observations. They are MBB mod¬ 
els with T = 40K, the DH02 template with a = 2.5, the CEOl 
template with L ~ 10®L©, and the SE_glx_l template of B13 
and C14 template. The minimum Lir required to observe them 
are explored by defining the complete Lir-z parameter space for 
the H-ATLAS, COSMOS, and GOODS-N surveys. The curve 
^completeis determined as follows: 

1. Calculate the minimum detectable IR luminosity at z for each 

SED (L[™(5(A),z)) from the SED template SiX) and the de- 


Table 4. The 3 ct depths of H-ATLAS, COSMOS, and GOODS, includ¬ 
ing confusion noise. 


Survey 

250 pm 

Depths[mJy] 
350 pm 

500 pm 

reference 

H-ATLAS 

18.0 

21.0 

27.0 

Smith et al. 

2011; 

COSMOS 

15.2 

17.0 

19.2 

Oliver et al. 

2012 

GOODS 

5.7 

7.2 

9.0 

Elbaz et al. 



tection limit (Tablej^ for each survey. The limit for each 
SED template is shown as colored curves in Eigure[^ 

2. Take the maximum of L[™(5 (A),z) at each z. The result 
^complete is the luminosity required to detect a source hav¬ 
ing any of the SEDs listed in Table in the three SPIRE 
bands (black curve in Eigure]^. 

The region above the black curve in each panel shows the 
Lir-z range where sources with these five SEDs are all detectable 
(complete for all SED shapes). 

We find that 27 of the 57 galaxies in our sample are within the 
H-ATLAS completeness limit. Eor the GOODS detection limits, 
the number of such galaxies increases to 52. In the worst case of 
the H-ATLAS detection limits, we have 19 galaxies at z ~ 3 and 
5 galaxies at z ~ 4 within the completeness limits, whereas z 5 
and 6 samples are not complete. Eor the deepest field, the non¬ 
detection of galaxies with these templates intrinsically indicates 
there are fewer or no objects with such SED shapes rather than 
the selection effect at z ~ 3 or 4, implying that the dust temper¬ 
ature is intrinsically high for high-z galaxies. In fact, if a galaxy 
with T — 30K is above the GOODS detection limit, the dust mass 
required is larger than 10*®M© at z > 4. It is not surprising that 
these kinds of galaxies are rare. 

Another question is whether galaxies that are fainter than the 
completeness limits are colder. Examining the properties of our 
UV-selected sample, we found that the EIR faint objects at z ~ 
4 are not located in the upper left region, suggesting that they on 
average have high dust temperature. Therefore, the UV-selected 
galaxies with cold dust temperatures may exist, but should be 
much rarer than galaxies with high dust temperatures. One might 
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Fig. 7. Comparison of SPIRE colors for UV-selected galaxies from Hei- 
|nis et al.| ( |2014[ l (blue) and our high-z sample (red) at z ~ 3 (upper panel) 
and z ~ 4 (lower panel). In each panel, the blue dots show the colors de¬ 
rived from the stacked UV sample from |Heinis et ah] l |2014^ .The red 
triangle indicates the average colors of our high-z sample. The error 
bar indicates the standard deviation. Filled black squares present colors 
for the galaxies from |Saintonge et ar]j2013) . The asterisks indicate the 
colors given by MBB models with dust temperatures of 30K, 45K, and 
60K. 


argue that galaxies may exist with no detectable UV nor FIR, 
however, the lensed galaxies in our sample are intrinsically faint, 
and there is no galaxy located in the upper left region of the 
color-color diagram. Therefore, it is likely that the cold dust 
galaxies are intrinsically few at high-z. The high dust temper¬ 
ature of high-z galaxies may be due to the strong star formation 
activity in these galaxies. 


However, the argument is not conclusive because our sample 
still has some potential selection biases toward luminous galax¬ 
ies. First, our requirements that the galaxies have spectroscopic 
redshifts can bias the sample to be brighter than the formal de¬ 
tection limits. This effect may not be significant since similar 
results can be found for the sample of |Swinbank et ah ( |2014| l, 
who used the ALMA LESS sample with photometric redshifts 
2 < z < 5; few galaxies are found with Ssoq/S^sq > 1.5 (see 
their Figure 5). Second, although we can explore galaxies fainter 
than the detection limits taking advantage of lensing magnifica- 
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H-ATLAS 



z 


Fig. 8. Completeness limits Ljr for detecting all SED shapes. The up¬ 
per, middle, and bottom panels show the complete limits for H-ATLAS, 
COSMOS, and GOODS-N surveys, respectively. In each panel, the 
galaxies of our high-z sample are plotted as diamonds (Lir from Ta¬ 
ble [TJ, and the Ljr for each SED template is plotted in colored line 
(red: B13; yellow: DH02 with a = 2.5; green: C14; blue: CEOl with 
Lir = 10®Lq; magenta: MBB with T = 40A'). The black line is the 
maximum of the colored lines, above which the observation should not 
be biased against these SEDs. 
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Fig. 9. Color-color diagram of the 70,000 mock galaxies generated from 
M12 templates. The color indicates the redshift of the mock galaxy. 
The dotted line shows the color cut given by |Amblard et al.] ( |2010j l. 
Galaxies above this line are expected to have z ~ 2.6 and those under it 
are expected to have z ~ 1.8. The dashed line shows the selection of 500 
^m risers. The upper left region separated by this line selects galaxies 
with z > 4. The solid lines are our color cut method according to the 
distribution of the mock galaxies. The solid lines divided the diagram 
into four parts. From lower right to upper left these parts select galaxies 
at z < 0.5, 0.5 < z < 1.5, 1.5 < z < 2.5 and z > 2.5. 


tion, the faintest galaxies in our sample have Ljr ^ 10'*L©. It is 
possible that cold dust exists in galaxies fainter than ^ IO^'L©. 


5.2. Sampling high-z galaxies using the SPIRE color-color 
diagram 


At high redshifts, since the available spectroscopic information 
is extremely limited, it would be useful to be able to select high-z 
candidates from photometric data. The SPIRE colors have been 
used for selecting high-z candidates in several studies. For exam¬ 
ple, some studies select 500 pm risers with ^500 > ^350 > ^250 
(e.g., |Riechers et al.|2013[ [Dowell et al.|2014[ ) as high-z candi¬ 
dates. Galaxies with ^500 > ^350 > ^250 are likely to be at red- 
shift z > 4 according to the shape of the IR SED. The galaxy 
HFLS3 at z ~ 6 was found through this method. Another method 
is to model the redshift distributions on the SPIRE color-color 
diagram and compare with observations (Amblard et al.|[ 20 T 0 


Roseboom et al. 


|2012|l. From the redshift distributions, |Am^ 


blard et al. (2010 1 also dehne a color cut with 8500/8250 > 0.75 


and < 0.75, which divides their H-ATLAS galaxy sample into 
z ~ 2.6 and z 1.8 populations. 


However, all these methods are based on MBB models, 
where the redshift evolution of the SEDs is not included. There¬ 
fore, more sophisticated SED templates should be used to inves¬ 
tigate the correlation between redshift and colors. We showed 
in Section 1^ that the Ml2 templates can reproduce the SPIRE 
colors that are consistent with the high-z observations. The M12 
library also gives the evolved SEDs for z < 3. Therefore, in this 
section we use the M12 library to re-investigate these methods. 


Using the MS templates of M12, we construct 70,000 mock 
galaxies withz G [0.0,7.0]. We added an extra Gaussian standard 
deviation of 10 % to the fluxes and then calculate the colors of 
these mock galaxies. In previous studies, 10% is commonly used 


Sample with photo-z 




z z 




Fig. 10. Redshift distributions of the four subsamples of the 783 galax¬ 
ies divided using the color-cut method defined in this work (Equations 


|2]to0. 


(e.g., Amblard et al. 2010 Robson et al. 2014| l. Using reasonably 
larger values ( 20 %) does not affect our conclusion. 

We overplotted the selection methods of 500 pm riser 
(dashed line) and [Amblard et al.] (2010; dotted line) in Figure]^ 
Both methods are based on MBB models. It seems that the 500 
pm riser method selects a very pure high-z sample. The color 
cut dehned by Amblard et al. (2010 1 is able to divide the mock 
galaxy sample into high-z (z ~ 2 . 6 ) and low-z (z 1 . 8 ) popu¬ 
lations, whereas the redshift range in each population is quite 
large. 

As shown in Figure [^z > 3 galaxies are located in a smaller 
area on the color-color diagram predicted by the M12 library 
comparing with that predicted by MBB models (see Figure 1 in 
[Amblard et al.|20l0| l. Based on the distribution predicted by the 
M12 library, we redehne the color cuts as follows: 



^<0.4 ^-0.3. 

‘^350 ‘^350 

( 2 ) 

‘^250 

- 0.3 < — < 0.65 — - 0.25 

(3) 

‘^350 

‘^350 ‘^350 

‘^250 

-0.25 < — < 1.1 — -0.15 

(4) 

‘^350 

‘^350 ‘^350 


‘^350 ‘^350 

(5) 


Equations]^ to 1^ dehne four regions R0(z < 0.5), Rl(0.5 < z < 
1.5), R2(1.5 < z < 2.5) and R3(z > 2.5) in Figure]^ separated 
by the three solid lines. 

In the case of our spectroscopically conhrmed high-z sam¬ 
ple (57 galaxies at z > 2.5), 40 galaxies are in R3 region. [Am- 1 
[blard et ^ ( 2010| l’s high-z color cut selects 42 of these galax¬ 
ies. The method of 500 pm risers selects 23 of these galax¬ 
ies. The completeness of the 500 pm riser selection is 40%, 
worse than both our method work and [Amblard et al.j ( |2010| l. 
Since the 500 pm riser selection method is designed for select¬ 
ing z > 4 galaxies, we also check the selection completeness for 
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SOOpm risers selection Amblard+2010 selection 




Fig. 11. Photometric redshift probability distributions of a sample of 
250 pm selected galaxies. Left panel: The filled histogram shows the 
redshift distribution for the subsample of 500 pm risers (5500 > 5350 > 
S 250 ) and the open histograms show the redshift distribution for the re¬ 
maining galaxies in the sample. Right panel: The filled histograms show 
the redshift distribution of the high-z subsample (z ~ 2.6) selected using 


Amblard et al. 1 2010|)’s color cut (S 500/5250 > 0.75), and the open his- 
tograms show the redshift distribution of the low-z subsample (z ~ 1.8) 
selected using 5500/^250 < 0.75. 


Table 5. Statistics of the redshift distribution (photometric z) in each 
region on the color-color diagram divided using Equations to and 
other methods. 



N 

Mean 

Median 

STD 

contamination 

All 

783 

1.31 

1.06 

1.00 


RO (z < 0.5) 

15 

0.18 

0.16 

0.09 

0% 

R1 (0.5 <z< 1.5) 

90 

1.09 

0.96 

0.94 

62% 

R2(1.5<z<2.5) 

463 

1.35 

1.08 

0.99 

83% 

R3 (z > 2.5) 

215 

1.39 

1.13 

1.03 

80% 

Amblard+(z > 2) 

384 

1.39 

1.12 

1.02 

77% 

Amblard+(z < 2) 

399 

1.23 

1.02 

0.98 

18% 

500 pm risers (z > 4) 

22 

1.13 

0.96 

1.09 

95% 


icantly lower than the expectation for this method. The contam¬ 
ination becomes higher than 50% at z > 2 for all the methods, 
and the dispersion of redshifts in all the z > 0.5 subsamples in 
Table|^is quite large (~ 1). 

The contamination is partly because of the sample selection. 
The requirement of association with 24pm source and an en¬ 
try in the SWIRE Photometric Redshift Catalog biases the sam¬ 
ple against sources with ^ 500 /S 24 > 200 as well as sources with 
z > 1.5 ( |Rowan-Robinson et al.|2008]|2014| l. Our requirement of 
detection at 250/4m further biases this sample to lower redshifts. 
The lack of z > 1.5 sources in this sample therefore causes poor 
statistics at z > 1.5. However, assuming galaxies without 24pm 
or 250/4m detection are all perfectly divided by the color-cut 
methods, we can estimate a minimum contamination rate, which 
is 53% for R3 (z > 2.5 cut), 53% for [Amblard et al.] ( |2010[ )’s 
z > 2 cut, and 50% for SQQpm risers cut. The minimum con¬ 
tamination is still high for high-redshift color cut. The error of 
photometric redshift and the cross-matching accuracy is also rel¬ 
evant to the contamination. The photometric redshift error fo r 
this sample is about 4% in (1 -fz) (Rowan-Robinson et al. 2014 1 , 
which increases with the redshift. At higher redshift, the associ¬ 
ation of submillimeter sources with 24pm sources also becomes 
less reliable. Apart from the above causes, there is also intrinsic 
dispersion caused by the variation of galaxy properties because 
we are considering an average SED in a given redshift bin when 
designing the method. 

From the discussion above, it seems that although the three 
methods can select samples with different mean redshifts, the re¬ 
sulting samples have either a large dispersion in z or miss a con¬ 
siderable number of galaxies 33%) at high redshift. There¬ 
fore, it is difficult to properly sample galaxy populations at dif¬ 
ferent redshifts using only the three SPIRE bands. 


z > 4 galaxies in our sample and found the completeness is only 
approximately 66% (14 of 21). 

To further test the three methods, we use the 500 /rm se¬ 
lected sample with photometric redshift information injRowam] 


Robinson et al. (2014; 783 galaxies, see Section 2.2 for details 


of this sample). The large size of this sample is necessary to test 
the color-cut methods because these methods are designed for 
statistical use. Figure shows the redshift distributions of the 
four subsamples of the 783 galaxies divided using the color-cut 
method defined in this work. The redshift distributions of sub¬ 


samples divided using the methods of 500 /4m risers and Am- 
blard et aL|p010|l are shown in Figure[^ The mean median red- 


shift, the standard deviation, and the contamination from other 
redshifts for each subsample are listed in Table 

From the statistics, we find that at z < 1.5 the color cuts de¬ 
fined in this work can separate galaxies into subsamples with 
different mean redshifts, and the Student’s t-test shows that the 
difference between these subsamples is significant (the index 
a < 0.05). The mean redshift of galaxies is 0.18 in RO region, 
and 0.96 in R1 region, consistent with our definition of these 
regions using mock galaxies. However, at z > 1.5 the mean red¬ 
shifts in R2 and R3 regions are lower than expected. The high- 
and low-redshift subsamples given by the method of j Amblard] 
et al. (2010i have mean redshifts of 1.43 and 1.23, respectively. 
The Student’s t-test shows that the difference is also significant. 
However, the mean redshifts are different from [Amblard et al.j 
( 2010| l values for H-ATLAS galaxies (2.6 for the high-z popula¬ 
tion and 1.8 for the low-z population). The method of 500 pm 
risers selects a high-z sample with mean redshift of 1.53, signif- 
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6. Conclusion 

We collected a sample of 57 galaxies with reliable redshifts and 
SPIRE flux measurements at z > 2.5, and compare the SPIRE 
colors of this sample with those derived from different SED tem¬ 
plates at different redshifts. The templates are taken from the 
pre-Herschel and post-Herschel SED libraries, including local 
templates of CEOl, DH02, and C14, and high-z templates of 
E11,M12, andB13. 

From the SPIRE color-color diagram, we find that ~ 95% 
of the galaxies at z > 2.5 have a SPIRE color 8500 / 80,50 < 1-5, 
implicating that the dust temperatures in these galaxies are high. 
Therefore, local SED templates may not be suitable for these 
galaxies due to dust temperature evolution. From the compari¬ 
son of the high-z data with the local templates, we find that the 
SPIRE colors given by the C14 template, which is based on lo¬ 
cal galaxies that are too cold to describe the average colors of 
our high-z sample. The local calibrations for CEOl and DH02 
libraries are also not able to describe the average SPIRE colors 
for high-z galaxies. However, both the CEOl and DH02 tem¬ 
plates can fit the SPIRE colors better when the parameter Lir of 
the CEOl templates is adjusted to ~ IO^'L© or the parameter a 
of the DH02 templates fixed to 1.5. Among the high-z templates, 
the templates of Ml2 give the most satisfying fit, implying that 
their assumption of the SED evolution is suitable for the high 
redshift. 

Comparing with UV-selected samples at z = 3 and 4, we find 
that the UV-selected sample have even higher dust temperatures 
than our sample. We discussed the possible bias of the samples 
due to the detection limits, and find that the higher than local 
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dust temperature for high-z galaxies is unlikely to be due to ob¬ 
servational biases. 

The efficiency of selecting high-z galaxies using the SPIRE 
color-color diagram is also discussed. We defined color cuts ac¬ 
cording to the MS template of M12 to sample galaxies at differ¬ 
ent redshifts, and compare the method with the method using 500 
/rm risers and the color cut defined by Amblard et al. ( |2010| l. We 
have shown that these methods can divide our 250 fim selected 
sample into subsamples with different mean redshifts., however, 
the dispersion in each subsample is quite large. Additional infor¬ 
mation is needed for better sampling. 
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Appendix A: SED tracks on SPIRE color-color 
diagram predicted by different templates 

Figure |A.1| plots color-color relations for the templates with 
Lir ~ 1W10^ 10'“, 10", 10'2 , and lO'^ L©. Ljr = lO"/.© 
templates, and Ml2 MS templates produce very similar SPIRE 
color-color relations. 

Figure |A.2| shows that the a = 2.0 template agrees well 
with M12 templates at the local universe. At z = 2, a must 
reduce to ~ 1.5 to agree with the evolution trend suggested by 
M12 templates. The evolution of M12 MS templates is similar 
to a variation of the value of a between 1.5 to 2.0. 


Figure A.3 shows the four templates selected from the 
B13 libraries. The template SF_glx_2 (the template for star¬ 
forming galaxies derived from COSMOS field) reproduces the 
SPIRE colors very similar to the M12 MS templates at z > 
4. The SE_glx_l template (star-forming galaxies derived from 
GOODS-S field) reproduces much colder colors than the M12 
templates. The MIRex_SE_glx (MIR excess star-forming galax¬ 
ies) and Obs_SE_glx (obscured star-forming galaxies) templates 
have hotter SPIRE colors than the M12 MS templates. 
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Fig. A.l. SPIRE colors of the templates of CEOl (black dotted lines with vertical bars) with Lir ~ 10*, 10®, 10*®, 10* ’, 10*^ , and 10** Lq . Red 
solid line with squares are the colors given by the MS templates of M12. 
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Fig. A.2. SPIRE colors of the templates of DH02 templates (black dotted lines with vertical bars) with a = 1, 1.5, 2, and 2.5. Red solid line with 
squares are the colors given by the MS templates of Ml2. 
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Fig. A.3. SPIRE colors of the selected templates of B 13 templates (black dotted lines with vertical bars). Red solid line with squares are the colors 
given by the MS templates of Ml2. 
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